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The one-dimensional Hubbard model is integrable in the sense that it has an infinite family of 
conserved currents. We explicitly construct a ladder operator which can be used to iteratively gen- 
erate all of the conserved current operators. This construction is difi^erent from that used for Lorentz 
invariant systems such as the Heisenberg model. The Hubbard model is not Lorentz invariant, due 
to the separation of spin and charge excitations. The ladder operator is obtained by a very general 
formalism which is applicable to any model that can be derived from a solution of the Yang-Baxter 
equation. 

PACS numbers: 71.10.Fd, 75.10.Jm 
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The one-dimensional Hubbard model has attracted 
considerable interest because it is one of the few examples 
of a model for strongly correlated electrons that is exactly 
soluble [0 . The fact that it describes a doped Mott insu- 
lator and exhibits spin-charge separation (i.e., spin and 
charge excitations are independent of one another) has 
been argued to be relevant to understanding the unusual 
metallic properties of high-temperature superconductors 
1^. Furthermore, the model has rich mathematical struc- 
ture it is solvable by both the co-ordinate Bethe 
ansatz Q and the algebraic Bethe ansatz it has a 
hidden so(4) symmetry and is integrable in the sense 
that it has an infinite family of conserved currents [Q. 
The latter is a consequence of the fact that the model 
can be derived from a solution of the Yang-Baxter equa- 
tion Integrability is not just of mathematical inter- 
est because it may have implications for dissipationless 
transport the coherence of interchain transport [p^ , 
and whether the energy level spacing follows a Poisson 
distribution or the Gaussian orthogonal ensemble distri- 
bution characteristic of quantum chaotic systems [ pT| . 
Furthermore, integrability has been essential to recent 
exact calculations of transport properties in mesoscopic 
electronic devices [ p^ . 

In this Letter we clarify the algebraic structure under- 
lying the integrability of the Hubbard model by using 
the Yang-Baxter equation to explicitly construct a sin- 
gle operator B (known as the ladder operator) which can 
be used in a simple recursion relation (equation (^ be- 
low) to generate the whole family {t''"''}^o conserved 
current operators, i.e., operators that commute with the 
Hamiltonian and one another. This result is surprising 
in light of the lack of invariance in the model under the 
lattice version of the Poincare group. 

For continuum field theories in 1-1-1 dimension the gen- 
erators of the Poincare group are B. P and H , being the 
generators of Lorentz boosts and translations in space 
and time, respectively. {P and H are also the total mo- 
mentum operator and Hamiltonian, respectively). They 



obey the closed algebra 

[B,H] = P, [B,P]^H, [H,P] = Q. (1) 

It is extraordinary that a wide range of integrable lat- 
tice models (including the Heisenberg | |l3|,p^ , Calogero, 
Toda, |l^ and supersymmetric t — J [[15[ models) are 
invariant under a generalization of the Poincare group 
involving the entire infinite set {t*-"-'}^o of conserved 
currents. They satisfy the algebra 
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where t^"-* and t'^' are the momentum operator P and 
Hamiltonian H, respectively. Here the boost operator 
acts as a ladder operator on the infinite sequence of con- 
served operators. For the XX Z model the boost opera- 
tor can be identified with an algebraic element of a lat- 
tice Virasoro algebra |l^. This turns out to be of great 
practical significance because it permits the use of vertex 
operators for the determination of the energy spectrum 
and calculation of correlation functions JTzt . 

A crucial property in the manifestation of Lorentz in- 
variance in the above models is the fact that the R matrix 
which is a solution of the Yang-Baxter equation (equation 
(^) below) has the difference property, R{u, v) — R{u — v) 
for the spectral parameters u and v. This is because the 
spectral parameter plays the role of rapidity variable. A 
uniform shift in both rapidity variables, corresponding 
to a change in the Lorentz frame, leaves the i?-matrix 
invariant. In this sense those solutions with the differ- 
ence property are invariant under a lattice version of the 
Poincare group ||l^ . 

In contrast, the Hubbard model is not Lorentz invari- 
ant ||l8| since it exhibits gapless excitations with differ- 
ent velocities. It is for this reason that spin and charge 
separate. As a result, the R matrix associated with the 
Hubbard model does not have the difference property and 
so its integrability is not as well understood. Although 
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Lieb and Wu Q] gave a co-ordinate Bethe ansatz solu- 
tion in 1968, it was not until 1986 that Shastry demon- 
strated the existence of an infinite family of conserved 
currents. This involved constructing a two-dimensional 
model in classical statistical mechanics with a transfer 
matrix that commuted with the Hubbard Hamiltonian 
0. It was achieved by mapping the model onto a pair 
of coupled spin chains using the Jordan- Wigner transfor- 
mation. Although Shastry conjectured that the i?-matrix 
satisfies the Yang-Baxter equation, it was some time be- 
fore a convincing proof was available Furthermore, 
the use of the algebraic Bethe ansatz method to repro- 
duce the solution of Lieb and Wu has only recently been 
achieved A significant consequence of this algebraic 
development is that it facilitates the use of the quantum 
transfer matrix method for the analysis of the thermody- 
namic properties at finite temperature PJl9[|. 

Grabowski and Mathieu claimed that there is no 
"matrix" ladder operator satisfying (^) for the model, 
motivating them to construct the first seven conserved 
currents by "brute force methods". We now show how 
for any model derived from a solution of the Yang-Baxter 
equation there is a one parameter family of ladder oper- 
ators B{v) such that the conserved currents satisfy (||). 
In cases where a solution to the Yang-Baxter equation 
has the difference property, the conserved currents have 
no dependence on the spectral parameter v. In this in- 
stance the construction for the ladder operator occurs as 
a particular case of the more general method we describe 
below. 

As an application of our general result we then consider 
a one-parameter of Hamiltonians which include the Hub- 
bard model as a special case (v = Q). Our approach has 
the further appeal that we work directly with the fermion 
operators of the model (rather than a two-dimensional 
statistical mechanics model) and that the so(4) invari- 
ance of the model is manifest throughout. 

The Yang-Baxter equation (or star-triangle relation) is 
central to exactly soluble models because it a sufficient 
condition for the validity of the Bethe ansatz [ pT|j22| . The 
corresponding equations for 1-1-1 dimensional quantum 
field theories are also known as the factorisation equa- 
tions because they imply that all possible decompositions 
of the N-particle scattering (S) matrix give the same re- 
sult as a product of two-particle S matrices . Consider 
a lattice model defined on L sites, each of which has a 
Hilbert space V . The matrix R{u, v) acts on the ten- 
sor product space V ^ V and satisfies the Yang-Baxter 
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i?i2(u, w)Ri3{u, v)R23{w, v) 

= -R23(w, v)Ri3{u, w)i?i2(u, w) 



(3) 



where the subscripts refer to the embedding of R{u, v) 
on the 3-fold space V (^V (SiV. From a solution to the 



this equation we define a transfer matrix 

t{u,v) = tro {RoLiu,v)...RQ2{u,v)Roi{u,v)) 

where tr denotes the trace over V (when is a superspace 
we use the supertrace, i.e., the trace over the bosonic 
states minus the trace over the fermionic states). It fol- 
lows from the Yang-Baxter equation (||) that 



[t(u, w), t{w, v)] — 
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for all values of the parameters u and w. 

An assumed feature of the R matrix is the regularity 
property, i.e., R{u,u) = P, with P being the permuta- 
tion operator. (A phase of (-1) is gained whenever two 
fermionic states are interchanged.) Using this property, 
the Hamiltonian is defined to be 



H{v) 



_i dt{u,v) 



du 
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where T = t(u,u) = Pil---PizPi2 is the translation op- 
erator. This yields 



L 

with the local Hamiltonian given by 

dR{u, v) 



(6) 



h{v) 



-P 



du 



Above and throughout periodic boundary conditions are 
imposed. For later use, it is convenient to consider the 
series expansion for the i?-matrix 

i?(w, v) = P{I +{v- u)h{v) + l/2{v - uff{v) 

+ l/(S{v-ufg{v) + ...) (7) 

Expressing the logarithm of the transfer matrix in a 
power series expansion 

lnt(u,z;) = ^i^i^i(")(«), (8) 

n=0 

it is apparent, in view of (|^), that 

[t("'(w), t("'(i;)] = 0, Vm,n 
and in particular 

Consequently, the Hamiltonian H{v) is integrable since 
the set of operators {t'^^Hv)} provide a set of conserva- 
tion laws for the system. Note that in the case where the 
i?-matrix does not have the difference property, we have 
the generic feature that the Hamiltonian and higher con- 
served charges will always have non-trivial dependence 
on the variable v, as can be seen from (p|). 
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We now construct the parameter-dependent ladder op- 
erator B(v). Differentiating the Yang-Baxter equation 
with respect to w, then setting w = v and premul- 
tiplying by the permutation operator Pj(j+i) yields an 
analogue of the Sutherland equation j2^] 

\^]U+i) ^o(j+i) (w, v)Roj{u, v)] 
= Rou+i){u,v) — — Roj{u,v). (9) 

An immediate consequence of (|) is that 

[Hiv), tiu,v)]=0 

which also follows from (||). 

The ladder operator B{v) is defined in terms of the 
local Hamiltonians hj = (w) through the relation 



(10) 



where j are the elements of the integers modulo L. A 
consequence of the generalized Sutherland relation (||) is 
now 



= 5Z hj+i + hj-i] - H' 
j 

3 

where the prime denotes a derivative with respect to v. 
The computation of t*-^' can be simplified by invoking the 
generalized Reshetikhin condition (cf. ||2^ ) 

[hi2 + h23, [hi2, ft.23]] 

= 2:23-2^12 (12) 

with the two site operator given by 

x = 2h^ + g-3hf + 2[h', h] - h" . 

The Reshetikhin relation is obtained by applying 
d'^/du^dw to (||) and using (^. Omitting the details, 
this yields the result 



(11) 



The relation (|11|) permits us to deduce the recurrence 
relation @) from (|). 

The definition of the ladder operator here is different 
from the cases considered in p"^-[T^] by the inclusion of 
the differential operator. This term is not required for 
those cases with the difference property since it is ap- 
parent from (||) that the conserved currents t*^"^ have no 
dependence on v and hence (^) still holds. However, in 
this instance (0) becomes 



[B, t{u - v)] 
which can be integrated to 



dt{u — v) 
du 



t{u + A) = exp(AB)t(w) exp(-AB). 

The parameter u characterizes the Lorentz frame for the 
transfer matrix t{u), showing that B is the generator of 
Lorentz boosts in this context . This is clearly not the 
case for the Hubbard model where Lorentz invariance is 
not present. 

From (^ the operators t^^\v) may be calculated iter- 
atively. We find the following expressions for the leading 
terms 

t(o) = inr. 
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+ Y(f'3~93~h,.h'^~2h'^.h,). 
3 

Above, c is a constant determined by the normalization 
of R{u,v). It can always be chosen to be zero. Before 
turning to the particular case of the Hubbard model we 
stress that the above construction of the ladder operator 
is valid for any solution of the Yang-Baxter equation (|^). 

Consider the four dimensional local Hilbert space V 
spanned by the states 

|0), IT), li), ITi). 

Introduce the spaces W{a) with basis {|0} , \a)}, a —1 1 
so that V = W{'\) X VFd). For each tensor space 
W{a)®W{a) there is a solution of the Yang-Baxter equa- 
tion (with difference property) given by 

TZijiu - w) = cos(w - w)(l - rijo- - rija) 

+ sm{u - v){nicr + rijcr - 2n,a-njcr) 



'4a^3<^ 



Here cj^. and Cja are the creation and annihilation oper- 
ators with spin u{=] , |) at site j and rij^ = '^ja'^ja the 
density operator. The associated Hamiltonian obtained 
through (^ is that for free fermions. 

It has recently been shown that the following R- 
matrix is also a solution of the Yang-Baxter equation 
acting on V 
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Ri] (u, v) 



Kin 



cos(u — v) 



tanh(e'(u) - e{v)) 



cos(m + v) 

X 7^^(^i + v)nl^{u + v){l- 2n,t)(l - 2n,x) (13) 
with 9(u) defined through the relation 

sinh26{u) = ^ sin 2m. 

An important consequence of (^) is that it aUows for 
the construction of a generalized Hubbard model (with 
spectral parameter dependence) as noted in ||2^. The 
identification of this generalized model is paramount in 
the construction of the laader operator. Explicitly, the 
local Hamiltonians read 

J 



E 



U 



4cosh26'(i;) 



(14) 



where 



- i3<^ 



(w) = cos^ u(l — 2nicr) — sin^ v{l — 2nj„) 



+ sin 2v{cl^Cj^ - cj^c,^). (15) 

It is clear that (|^) with (|l^) reduces to the usual Hubbard 
model when v — 0. In particular, we find 

/i:,.(0) = U/2 [(1 - 2n.T)(4^c,i - c^c^i) 
+ (1 - 2n,i){cl^cj^ 



C]|C,|) 



(16) 



which plays an important role in the explicit construction 
of the higher conserved operators discussed previously. It 
is important to note that there is no way of determining 
h'{0) directly from the usual Hubbard model. 

Substituting g (which is obtained through (Q) and (|li 
and eqs. (|5 16), all evaluated at u = 0, into the above 
expressions for t^^' and t^^\ we recover the expressions 
found previously for the first and second |^,^ non- 
trivial conserved currents (modulo a constant term and 
multiple of H). A well known feature of the Hubbard 
model is that all the integrals of motion, except for the 
translation operator, are invariant with respect to the 
so(4) Lie algebra when the lattice length is even [ p5| . 
Significantly, the ladder operator (|lO|) in the case of the 
Hubbard model is also so(4) invariant for an even length 
lattice. 

To conclude, a systematic method for obtaining the 
conserved currents in the Hubbard model has been de- 
scribed which employs the use of a ladder operator. We 
emphasize, however, that the construction presented here 
is entirely general and may be applied to any Yang- 
Baxter integrable system. 
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